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Abstract 

We demonstrate the existence of a simple physical picture of superconductivity for extremely over- 
doped Cu02 planes. It possesses all characteristic features of HTS, such as a high superconducting 
transition temperature, the d x 2_ y 2 symmetry of order parameter, and the coexistence of a single elec- 
tron Fermi surface and a pseudogap in the normal state. Values of pseudogap are calculated for different 
doping levels. An orbital paramagnetism of preformed pairs is predicted. 

PACS numbers: 74.20.-z; 74.78.-w 

1 Introduction 

In this work (see also earlier Letter [1]), we demonstrate that, in the phase diagram of cuprate HTS, a small 
region exists in which the characteristic features of HTS can be easily understood on the bases of a simple 
theory. The above characteristic features are such as a high superconducting transition temperature, the 
d x 2_ y 2 symmetry of order parameter (see [2]), and the coexistence in the normal state of a single electron 
Fermi surface and a pseudogap [3] . The latter phenomenon is usually attributed to the presence of preformed 
(i.e., in the normal state) electron pairs (in particular, bipolarons [4-8]). 

The aforementioned small region in the phase diagram is situated in the vicinity of the maximal hole- 
doping level x = x c compatible with superconductivity. The superconducting transition temperature T c is 
zero for x > x c , so it is low in our region near x = x c . However, T c increases with decreasing x for x < x c 
in such a way that at the boundary of the region (i.e., for x c — x ~ 1) it is quite high. 

Two features of our small region are important to make a simple physical picture possible. These are 
relatively low T c and the clear nature of the normal state as mostly the conventional Fermi liquid. 

We calculate the pseudogap. With increasing x, the pseudogap decreases for x < x c . As well as T c , the 
pseudogap disappears at x = x c . However, for larger doping levels x > x c , it appears again. 

As a new prediction, we show the existence of an unusual orbital paramagnetism of the preformed (singlet) 
pairs, which probably can be experimentally separated from the Pauli spin paramagnetism of single electrons 
and the Landau diamagnetism of single electrons and pairs. 

2 Pair quasiparticles 

The key point is the existence of very mobile pair quasiparticles in crystals in conditions of the tight-binding, 
i.e., if the energy of electron-electron interaction at a distance on the order of atomic spacing considerably 
exceeds the electron-tunneling amplitude to neighboring lattice cites. Quasiparticles of this type were studied 
earlier [9] in helium quantum crystals and more recently by Alexandrov and Kornilovitch [7] as a model of 
bipolarons in HTS (see also [10]). 

Let us consider two electrons localized at neighboring (1 and 2 in the figure) copper atoms (to be more 
precise, in unit cells containing these atoms) forming a square lattice in the Cu02 plane. The electron 
tunneling from 2 to 4 or 6 does not change the energy of the system in view of the crystal lattice symmetry. 
The same is true for the electron tunneling from 1 to 3 or 5. Owing to this type of transitions, an electron 
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Figure 1: C11O2 plane: (•) Cu atoms and (x) O atoms. 

pair can move as a whole over the entire plane, since the 2 — > 4 transition can be followed by the transition 
1 — ► 7 or 1 — ► 3, and so on. Since the transitions do not change the energy of the system, the motion is fully 
coherent. An electron pair behaves as a delocalized Bosc quasiparticle. 

To calculate the quasiparticle spectrum, we consider the localized states of a pair, 

\r,r',ai3)=ct a c+ \O), (1) 

where c+ Q are the electron creation operators at point r with spin projection a =t , | and |0) is the electron 
vacuum. 

The effective tunneling Hamiltonian H e f f is defined by the matrix elements of the operator 

H = tJ24 a c™, (2) 

rr'a 

which correspond to the transitions of one of the electrons to copper atoms that are next-to-nearest neighbors 
of the initial atom, in such a way that the energy of the system of two electrons remains unchanged. Here, 
t is the tunneling amplitude, which is known to be positive (see [2], p. 1004). 

Let a„ (n = x, y) be the square-lattice periods directed from point 1 to point 2 and from point 1 to point 
4, respectively. We have 

H eff |r,r + a x ,a/?) = t(\r + & x + a v , r + a x , a(3) + |r + a x - a y , r + a x , a/3) + 

|r,r + a^, a/3) + |r,r - a y ,a/3)) = t(- |r + a x , r + a x + a y ,(3a) + |r + a x - a y ,r + a x ,a/3) + 

|r,r + ay, a/3) - |r - a y ,r, (3a)), (3) 

where we used the antisymmetry of quantities (1) with respect to arguments (r, a) and (r', (3). Analogously, 

H eff |r,r + a y ,af3) = t(- |r + a y ,r + a x + a y ,[3a) + |r - a x + a y ,r + a y ,a(3) + 

\r,r + a x ,a[3) - |r - a x , r, /?«)). (4) 

The complete set of localized states of an electron pair is determined by the state vectors 

\r,n,a(3) = |r, r + a„, a/3) , (5) 

where r labels unit cells of the square lattice. 

The problem obviously splits into two independent problems for singlet and triplet pairs that are charac- 
terized by quantities (5), respectively, antisymmetric and symmetric about the spin indices a, (3. Assuming 
that the required stationary states of a pair are superpositions of localized states, 

]>>ge* kr | W /3) (6) 

r,n 
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with coefficients independent of r (this corresponds to a definite quasimomentum k), we obtain 

(.E(k) - e )V (x) = ± ± e*"*), 

(E(k) - e )V (!/) - t^ x \l ± e 4K *)(l ± e"^), (7) 

where the upper or lower sign corresponds to a singlet or triplet state, respectively. The conditions for the 
existence of a nontrivial solution ip( x \ tp^ to system (7) defines the energy E(k) of a pair quasiparticlc. 
Here eo is the energy of the initial localized state; n x = ka x and k v = ka^. Everywhere in formulas (7), wc 
omitted identical spin indices a(3. 

The minimal energy e m = min£'(k) = eo — At of a singlet pair is attained for k x = n y = 0. The same 
minimal energy of a triplet pair is attained for a nonzero quasimomentum k x = n y = it. This degeneracy is 
removed by taking into account the electron exchange in the initial localized pair. It is well known that this 
exchange is of an antiferromagnetic nature and, hence, singlet pairs possess the minimal energy. 

Thus, solitary Bose quasiparticles can exist in the Cu02 plane; these particles are characterized by a 
doubled electric charge and by zero momentum and spin in the ground state. It can readily be seen from Eqs. 
(7) that the effective mass of quasiparticles is m — h 2 /ta 2 , where a = |a x | = a y |. In addition, quasiparticles 
possess a specific quantum number n — x,y, which determines the orientation of a two-electron "dumbbell". 
Substituting E(k) = e m and k = into Eqs. (7), we obtain ip^ = —tp^ in the ground state. Since 
orientations n — x and n = y are transformed into each other upon the lattice rotation through an angle of 
it/2 and upon the reflection in the diagonal plane passing through points 1 and 3 in the figure, the ground- 
state wave function ip = ip^ — —tpiv) of quasiparticles transforms in accordance with the nontrivial ID 
representation (usually denoted by d x 2_ y 2) of the symmetry group of CuC>2 plane (see [2]). 

3 Superconductivity 

We further assume that all other two-, three-, etc., electron configurations localized at distances on the 
order of atomic spacing are energetically disadvantageous as compared to the pair configuration considered 
above. In addition, we assume that electrons are repulsed at large distances such that the electron-electron 
interaction energy is on the order of the one-electron tunneling amplitude. Under these conditions, only 
single electron Fermi particles and the pair Bose particles considered above play a significant role. 

Finally, let us assume that the minimal energy e m of pair quasiparticles is such that e m /2 is within 
a single electron energy band. We note the following. In conditions of the tight-binding, there are two 
different situations in which e m /2 can be within a single electron energy band. First, if single electrons and 
electrons in pairs correspond to the same energy band, single electron tunneling amplitude should be of the 
order of the electron-electron interaction energy in the pairs, while one-electron tunneling amplitude t in the 
pairs, introduced in Section 2, should be much smaller than the interaction. The last condition which is the 
condition of applicability of the procedure used in Section 2, can be result of the large polaron effect in the 
pairs. In the second situation, single electrons and electrons in pairs correspond to different bands. Both 
one-electron amplitudes can be of the same order in this case. The analysis carried out by Alexandrov and 
Kornilovitch in [7] shows that the conditions formulated above are likely to be realistic. 

Let us trace the change of the state of the system at T = upon an increase in the number of electrons 
(decrease in the hole-doping level). Until e m /2 > ep, only single electron quasiparticles arc present and the 
system behaves as an ordinary Fermi liquid. The condition e m /2 — ep determines the minimal hole-doping 
level compatible with the state of a normal Fermi liquid. Denote by n c the corresponding electron density 
n. Upon a further decrease in the hole-doping level, all additional n — n c electrons will pass into a Bose- 
Einstein (BE) condensate of pair quasiparticles (we everywhere consider the case of small n — n c values, 
for which the concentration of pairs is low and their interaction can be disregarded). The system becomes 
a superconductor. The superconducting order parameter represents the boson ground-state wave function 
ip = ^>( x ) normalized by the condition \ip\ 2 = (n — n c )/2; wave function transforms in accordance with the 
d x 2_ y 2 representation of the symmetry group of CuC>2 plane. 

It is important to note the following. In the system ground state (i.e., for complete filling of all fermion 
states with an energy smaller than ep), the uncertainty in the energy of a boson quasiparticle with low 
excitation energy e = k 2 /2m, arising due to its collisions with single electron Landau quasiparticles, is 
proportional to e 2 . As in the conventional theory of Fermi liquid, this is due, first, to the fact that the 
density of fermions in the vicinity on the order of e near ep with which the given boson can collide due to 
energy conservation, is low. Second, the statistical weight of the final states, to which fermionic transitions 
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are possible, is small. The probability of the boson decay into two fermions per unit time is also small: as 
suggested in the beginning of this Section, the boson must overcome a significant energy barrier. Thus, the 
proposed picture of superconductivity in the vicinity of maximal doping level remains valid even in the region 
of appreciable densities of fermions, where the interaction between bosons and fermions is significant. The 
critical electron density n c is determined from the condition that the electron chemical potential is equal to 
half of the minimal boson energy In the general case, the latter is a functional of the distribution function 
for single electron Landau quasiparticles. 

In calculating the superconducting transition temperature, the fermion distribution function may be 
considered as corresponding to T = 0, since the temperature corrections (proportional to T 2 ) to the ther- 
modynamic functions of Fermi liquid are considerably smaller than the corrections included below. 

The density of uncondensed bosons at a finite temperature T < T c is 

AT . f 2Mk 1 mT T 

The integral in Eq. (8) diverges at small e, so that it is cut off at e ~ r, where r is a small tunneling 
amplitude of electrons in the direction perpendicular to the CuC>2 plane. 

The excess number n — n c of electrons in the system is equal to the doubled sum of N' and number N 
of bosons in the condensate. This leads to the dependence of the superconducting transition temperature 
on the doping level for small values of n — n c 

mT c T c 

n-n c = —j log — (9) 

and the number of pairs in the condensate 

N °- — V-T e ]^)' (10) 

which determines the modulus of order parameter \ip\ 2 — Nq at finite temperatures. The superconducting 
transition temperature defined by Eq. (9) is quite high. To within the logarithmic term, this temperature 
is on the order of one-electron tunneling amplitude t at the boundary of the applicability region (i.e., for 
n — n c ~ a~ 2 ). The possibility that the superconducting transition temperature may have such an order of 
magnitude was pointed out in the aforementioned paper by Alexandrov and Kornilovitch [7]. 

The interaction of fermions with the BE condensate (effective electron-electron interaction) that is de- 
scribed by the order parameter tp creates an effective potential At acting on fermions as in conventional 
superconductors: 

H rnt = ^(A k c+ C ; ki + h.c). (11) 

k 

In view of the symmetry of ip, we have 

A k = V{k 2 x - k 2 y )^ (12) 

where k = k/|k| and V is invariant under the symmetry group. 

Owing to this interaction, fermions in the superconducting state considered acquire features typical of 
an ordinary superconductor with the d x i i symmetry. 



4 Normal state thermodynamics. Pseudogap 

The total number of pairs for T < T c is independent of temperature and equal to (n — n c )/2. The electron 
chemical potential for T < T c is also temperature independent and equal to fi — n(n c ) = e m (n c )/2 where 
e m = e m (n) is the pair minimal energy which depends on fermion density, as it is shown above. 

For T > T c , the fermion distribution function, as above, corresponds to T = 0, but with the temperature 
dependent chemical potential. The pair energy spectrum is E = e m (fi) + e where e = k 2 /2m. The pair 
density above T c is given by 

f°° 2-Kkdk 1 _ mT 1 
N ~J (2nh) 2 e^+O/T - 1 ~ l0g 1 - e -C/T ' ^) 
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The parameter £ (( 3> t) is defined by 

C=^>-2J/z (14) 
where Sfi = fx — /ti(n c ). With changing temperature, the total electron number conservation gives 

ra - n c = 2iV + — fyt. (15) 
o/j, 

From the last equation, we find ( = ((T) and then all other quantities. 

For n> n c and not too high temperature T -C T c log(T c /r), the pair density is determined by 

N(T) - N(T C ) dn/dfi 



N(T C ) 2(2 - de m /dfi) 

where N(T C ) = (n - n c )/2 and 



Te-^/ T (16) 



A p = T c log^ = — (n-n c ) (17) 
r m 

is the pseudogap for n > n c . As well as T c , it is zero at the critical value of the doping level n = n c . For 
higher doping level n < n c , (T c = 0), we have 

^■S^ (18) 

where 

4 -(*£-£)"--> < 19 > 

is the pseudogap for n < n c . The equation (18) takes place in the low temperature region T <C A p . For 
n < n c , the pseudogap A^ is the gap in the energy spectrum of the pair quasiparticles. For high temperatures 
T » A p , A' p (but Tct), the pair density is a linear function of temperature: 

where z is the solution of the equation Xz = e~ z with 

^—y^%r- (2D 

to 2 — oe m /Ofj, 
The entropy of pairs is determined by the equation 

to f 00 

where / = {e( e+ ^/ T — l} 1 . For n > n c in the low temperature region T -C A p , we have 

'e- A ^ T (23) 



(22) 



S(T) fS\ m A p ^_ Av/T 



T J t=t 2irh T 

where 

' S\ irm 



T J rp rp 12?l 



(24) 



T=T e 

S(T) is almost linear in T, with exponentially small deviations. For n < n c , the pair entropy is exponentially 
small at low temperatures T< AJ,: 

S (T) = ^L e -^' T . (25) 
At high temperatures T> A p ,Aj„ the entropy is 

= ^T. (26) 

The temperature independent factor a is determined by 

xdx 
e x - 1 

The entropy is again a linear function of temperature. 



Xz 2 . (27) 
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5 Orbital paramagnetism of pairs 

In this Section, wc show that the orbital motion of electrons inside the pairs cause a peculiar paramagnetism. 
Let a pair is at rest as a whole. For singlet pairs at k = 0, the Hamiltonian (3), (4) can be written as a 2x2 
matrix 

H = U (° I) =4fcri, (28) 



.1 

acting to the state vector i[>, 

^ = ^)(J)+^")(Jj, (29) 

where ip( n \ n = x,y, are quantum amplitudes of two orientations of the two-electron dumbbell, and o\ is a 
Pauli matrix. 

In the ir-state, coordinates of two electrons (with respect to the center of gravity of the pair) are x\ = 
—a/2, yi = and x 2 = a/2, y 2 = 0, respectively. In the y-state, we have x\ = 0, y\ = —a/2 and x 2 = 0, 
y2 = a/2. From this we find operators of coordinates for both electrons: 

xi = -x 2 = --\ n J, yi = -y 2 = --( . (30) 



2 \0 0/ ' yz 2 \0 1 

The operators of velocities are determined by the commutators 



n 

Simple calculation gives 



h,2 = UH,r h2 }. (31) 



2at fO -i\ 2at 
xi = -x 2 = -yi = Wi = — t- I ■ n = — ^ CT 2- (32) 



h \i ) h 
The operator of the pair magnetic moment which is directed along z-axis, is 

e x , . eta 2 

ij, = h z = — zJ^y - v x ) = — ^-^2 (33) 

where e is the electron charge and c is the velocity of light. 

In the presence of an external magnetic field B = B z , the Hamiltonian of the pair is 

H = Ata 1 - fj,B. (34) 

The eigenvalues of energy are 

E = e T 4t 



, («* " 211/2 



In weak fields, the minimal energy is 



2„4 



e a 



8h 2 c 2 

The average magnetic moment of the pair is 



where 



(35) 



E m m = eo - 4t — t—^—B . (36) 



W = -^P=aB, (37) 



r* = 7— : 5 ( 38 ) 



4fr 2 c 2 4mc 2 
is the pair paramagnetic polarizability. 

We note that pairs with k = in the upper energy band (lower sign in (35)) are diamagnetic. 
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The pair contribution to the paramagnetic susceptibility of 3D sample is 



X = 



e 2 a 2 
Amc 2 



N^) (39) 



where = N(T)/L is 3D density of pairs and N(T) is 2D density which is determined by formulas (16), 
(18), and (20). Here L is the distance between neighboring CuC>2 planes. 

Generally, we have three competing contributions to the magnetic susceptibility: orbital paramagnetism 
of pairs considered above, Pauli spin susceptibility of single electrons (pairs are singlet), and Landau dia- 
magnetism of single electrons and pairs. Spin susceptibility is isotropic. Orbital paramagnetism and Landau 
diamagnctism are both strongly anisotropic (magnetic moment is directed along z-axis independently of the 
direction of magnetic field) because of 2D character of single electrons and pairs. However, Landau dia- 
magnetism, especially in 2D case, is very sensitive to inhomogeneities. For example, it is easily suppressed 
by a localization of charge carriers. Orbital paramagnetism is finite at zero velocity of a pair as a whole. 
So, it has to be much more stable against inhomogeneities. We hope that orbital paramagnetism can be 
experimentally separated from two other contributions to the susceptibility. 
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